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Abstract 

We present a simple method to decompose the Green forms corresponding to a large 
class of interesting symmetric Dirichlct forms into integrals over symmetric positive semi- 
definite and finite range (properly supported) forms that are smoother than the original 
Green form. This result gives rise to multiscale decompositions of the associated Gaussian 
free fields into sums of independent smoother Gaussian fields with spatially localized 
correlations. Our method makes use of the finite propagation speed of the wave equation 
and Chebyshev polynomials. It improves several existing results and also gives simpler 
proofs. 

1 Introduction and main result 

1.1 The Newtonian potential 

Let us place the result of this paper into context through an example. Consider the Newtonian 
potential, the Green's function of the Laplace operator on R d given by 

fM-( rf -2) (d>3) 
$(x) = C d { 11 y - > for all x G R , x ^ 0. (1.1) 

\logl/|x| (d = 2) 

For d > 3 and any measurable function tp : R — >■ M such that t d ~ s ip(t) is integrable, the 
Newtonian potential can be written, up to a constant, as 

roc jj. 
\ x \-(d-2) = / t -{d-2) ^\ x \/ t ) — for all x£R d ,x^0. (1.2) 
Jo t 

This is true because both sides are radially symmetric and homogeneous of degree — (d — 2), 
where homogeneity of the right-hand side simply follows from the change of variables formula. 
In particular, ip can be chosen smooth with compact support and such that <£>(|cc|) is a positive 
semi-definite function on M. d . The last condition means that y(|x|) is positive as a quadratic 
form: for any / G C^°(M. d ) : that is, / : M. d — > R smooth with compact support, 

**(/,/):=/" <p(\x-y\/t)f(x)f(y)dxdy>0. (1.3) 
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Similarly, if d = 2, and (p : R — >■ R is any absolutely continuous function with <p(0) = 1 
and such that (f'(t) is integrable, then 

f°° dt 

logl/|x|=/ {<f{\x\/t) -<p(l/t)) — for all x £ R 2 , x / 0. (1.4) 
Jo * 

Indeed, for x ^ 0, 

logl/|x| = y>(0)logl/|x| =- / tp'(s)logl/\x\ds= iff {a) I -da, (1.5) 

jo io ■'■s/NI 

and thus, since is integrable, by Fubini's theorem, 

f°° ft\x\ dt poo dt 

iogi/N = y y ^^ ds 7 = J (v(*N) - f , (i-6) 

showing (|1 .4[) after the change of variables i i— >■ 1/t. Now suppose again that c/? is chosen 
such that y(|x|) is a positive semi-definite function on M 2 . Then the function M 2 3 x t— > 
(f(\x\/t) — (p(l/t) is positive as a quadratic form on the domain of smooth and compactly 
supported functions with vanishing integral: 



MfJ)--= (f(\x-y\/t)-tp(l/t))f(x)f(y)dxdy (1.7) 

<p(\x-y\/t) f{x)f(y) dxdy>0 



for all / G C^°(1R 2 ) with / / dx = 0. 

The above shows that the Newtonian potentials (jl.ip admit decompositions into integrals 
of compactly supported and positive semi-definite functions, with the appropriate restriction 
of the domain for d = 2. 

Let us only remark at this point that the positivity of a quadratic form has the important 
implication that it entails the existence of a corresponding Gaussian process, discussed briefly 
in Section Tl. 41 It is however also of interest in mathematical physics for different reasons [22J. 



1.2 Finite range decompositions of quadratic forms 

It is an open problem to characterize the class of positive quadratic forms, S : D (S) x D(S) — > 
R, that admit decompositions into integrals (or sums) of positive quadratic forms of finite 
range: for all f,g £ D(S), t > 0, 

( f°° dt 

S(f,g)= S t (f,g)-, 
Jo 1 

S t : D(S) x D(S) -> R, (1.8) 

W,/)>o, 

k S t (f,g) = if d( S upp(/),supp(g)) > 6(t), 

where 6 : (0, oo) — > (0, oo) is increasing and d is a distance function. The condition of finite 
range, the last condition in (jl.8p . generalizes the property of compact support of the function 
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(p in (jl.3p to quadratic forms that are not denned by a convolution kernel. The difficulty 
in decomposing quadratic forms in such a way is to achieve the two conditions of positivity 
and finite range simultaneously. Note that by splitting up the integral, one can obtain a 
decomposition into a sum from (jl.8p . and conversely, a decomposition into a sum can be 
written as an integral (without regularity in t). 

For applications, not only the existence, but also the regularity of the decomposition 
(|1.8p is important. Let (X,/j,) be a metric measure space, i.e., a locally compact complete 
separable metric space X with a Radon measure fi on X with full support (i.e., [i is strictly 
positive), C C (X) the space of continuous functions on X with compact support, and Cb(X) 
the space of bounded and continuous functions on X. Let us say that the decomposition (|1,8I) 
is regular if C c (X)nD(S) is 5-dense in D(S) and if every St has a bounded continuous kernel 
a t e C b (X x X): 

S t (f,g) = J St (x,y)f(x)g(y)dn(x)dfx(y) for all f,g€ C C (X) n D(S). (1.9) 

For the decompositions (jl.2p , (jl.4p , the kernels are of course given in terms of the smooth 
function ip by the explicit formula 

fa(x,y) =r^" 2) ip(\x-y\/t) for all x,y G R d , t > 0. (1.10) 

Note that for d = 2 the second term in (|1.4|) could be omitted by (|1.7|) . with the understand- 
ing that the quadratic form is restricted to functions with vanishing integral. It follows in 
particular that 

\M X >V)\ ^ Ci~ (d ~ 2) uniformly in all s,y G R d . (1.11) 

This reflects the decay of the Newtonian potential. Moreover, for all integers l x ,ly > 0, the 
derivatives of the kernel St decay according to 

\D l iD l ^ t {x,y)\ < dt-V-^r 1 *- 1 * , (1.12) 

reflecting that \D l $(x)\ < Q\x\~^ 2 ^ for all x G R d , x ^ 0. 

The main result of this paper is a rather simple construction of decompositions (|l,8p with 
estimates like (jl.lip for quadratic forms that arise by duality with Dirichlet forms in a large 
class. Let us call such forms Green forms motivated by the Newtonian potential, or Green's 
function, that is a special case. This is explained in Section 11.31 

The main idea of our method is that (I1.8P can be achieved by applying formulae like (jl.2p 
to the spectral representation of the Green form, and then exploiting finite propagation speed 
properties of appropriate wave flows. These are generalizations of the fact that if u(t,x) is a 
solution to 

d 2 u - Au = 0, u(0, x) = u (x), d t u(0, x) = (1.13) 
with compactly supported initial data uq that then 

supp(n(t, •)) C JVt(supp(«o)) (1.14) 

where N t (U) = {x G X : d(x, U) < t] for any U C X. 

The idea of exploiting properties of the wave equation in the context of probability theory 
is not new. For example, Varopoulos [3l] has used the finite propagation speed of the wave 



3 



equation to obtain Gaussian bounds on the heat kernel of Markov chains, by decomposing it 
into compactly supported pieces. Our objective is slightly different in that we are interested 
in the constraint of positive definite decompositions. 

Decompositions of singular functions into sums or integrals of smooth and compactly 
supported functions have a history in analysis. For example, Fefferman's celebrated proof of 
pointwise almost everywhere convergence of the Fourier series [T7j uses a decomposition of 
1/s on R like (jl.2p . albeit without using positive semi-definiteness. Hainzl and Seiringer |22j . 
motivated by applications to quantum mechanics such as [18] . decompose general radially 
symmetric functions, without assuming a priori that they are positive definite, into weighted 
integrals over tent functions. These, like y(|ic|) in (|1.2p . are positive semi-definite. They state 
sufficient conditions for the weight to be non- negative, and thus obtain decompositions like 
(11. 2p for a class of radially symmetric potentials including e~ m \ x \ on M 3 . Special cases and 
similar results have also appeared in earlier works of Polya [27] and of Gneiting [20, 21 . 

These results, like fll .2f) . make essential use of radial symmetry. One example of particular 
interest for probability theory — where radial symmetry is not given — is the Green's function 
of the discrete Laplace operator: 



Brydges, Guadagni, and Mitter [6] showed that also in this discrete case, the corresponding 
Green's function, or more generally the resolvent, admits a decomposition like (|1.8p into a sum 
(instead of an integral) of positive semi-definite lattice functions with estimates analogous to 
(|1.12p . Brydges and Talaczyck [10] gave a related construction which applies to quite general 
elliptic operators on domains in M. d , but estimates on the kernels of this decomposition are 
only known when the coefficients are constant. Their construction was adapted by Adams, 
Kotecky, and Miiller [1] to show that the Green's functions of constant coefficient discrete 
elliptic systems on Z d admit decompositions with estimates analogous to ()1.12p and that the 
decomposition obtained this way is analytic as a function of the (constant) coefficients. These 
results are all based on a constructions that average Poisson kernels. 

Our method, as briefly sketched earlier, is different from that of [Bl [T0| [U [l] and yields 
simpler proofs of their results about constant coefficient elliptic operators — both in discrete 
and continuous context. It furthermore naturally yields a decomposition into an integral 
instead of a sum (with integrand smooth in t), and gives effective estimates for decompositions 
of Green's functions of variable coefficient operators. 

1.3 Duality and spectral representation of the Green form 

Let us now introduce the general set-up in which our result is framed more precisely. For 
motivation, we first return to the quadratic forms defined by the Newtonian potentials f 1 1 . 1 [) : 




(1.15) 



e gZ d :|e|i=l 




(1.16) 



where 




(1.17) 
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These quadratic forms are not bounded on L 2 (R rf ), as is most apparent when d = 2. They 
are closely related to the Dirichlet forms given by 

E(u,v):=[ Vu-Vvdx, u,veC™(M d ). (1.18) 

JR d 

The correspondence between the two is duality: for all / G D(&), 

$(/,/)= sup jy /udz : u G C™(R d ),E(u,u) < lj . (1.19) 

This set-up admits the following natural generalization: Let (X, //) always be a metric 
measure space and L 2 (X) be the Hilbert space of equivalence classes of real- valued square fi- 
integrable functions on X with inner product (u, v) = (u, v)i?. Let E : D(E) x D{E) — > R be 
a closed positive quadratic form on L 2 (X) with D(E) C L 2 (X) a dense linear subspace. It is 
sometimes convenient to assume that E is regular, i.e., that C c (X)f]D(E) is F-dense in D{E). 
That E 1 is closed means that D(E) is a Hilbert space with inner product E(u,v) + m 2 (u,v)i2 
for any m 2 > 0. For the example (|1.18p . the domain of the form closure D(E) of C£°(M d ) is 
the usual Sobolev space H l (M. d ) and (u,v) H i = E(u,v) + (u,v) L 2 is the usual Sobolev inner 
product. 

It follows [29] from closedness that E is the quadratic form associated to a unique self- 
adjoint operator L : D(L) — > L 2 (X), 

E(u, v) = (u, Lv) for u G D(E), v G D(L), (1.20) 

where D(L) C D(E) is a dense linear subspace in L 2 (X). Moreover, self-adjointness of L 
gives rise to a spectral family and functional calculus. This means in particular that for any 
Borel measurable F : [0, oo) —> R, there is a self-adjoint operator, denoted F(L) : D(F(L)) — > 
L 2 (X), where 

POD 

F{L) := / F(A) dP A , (1-21) 
JO 

D(F(L)) := ju G L 2 (X) : F(A) 2 d{u,P\u) < ool (1.22) 

with P\ the spectral family associated to L, and (u, Pa^) is the spectral measure associated 
to L and n G L 2 (X). In these terms, E has the representation 

E(u,u) = \\L^u\\ L 2 {x) = / \d(u,P x u), u e D{E) = D(L^), (1.23) 

■/ spec(L) 

where E(u, v) for u 7^ u is defined by the polarization identity. Similarly, the corresponding 
Green form can be defined by polarization and 

$(/, /) = \\L-^f\\ L2{x) = I A" 1 d(u, P x u), f G £>(*) = D(X-I). (1.24) 

1/ spec(L) 

This representation will be our starting point for the decomposition of the Green form. Be- 
fore stating the result and its proof, let us sketch how the decomposition problem arises in 
probability theory. 
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1.4 Gaussian fields and statistical mechanics 

Even though the linear space D(E) is complete under the metric induced by the inner product 
E(u, v) + m 2 (u, v)i2 for any m 2 > 0, it is generally not complete for m 2 = 0. It may however 
be completed to a Hilbert space abstractly; we denote this Hilbert space by (He, (•, -)e)- 
Similarly, we can complete the domain D(&) to a Hilbert space under the quadratic form 
this Hilbert space is denoted by (H$, (■, •)$). He and H$ are dual in the following sense: The 
L 2 inner product can be restricted to 

(-,-): £>($) x D(E) -> R, (/, u) = (f, u) = (L~*f, l£u) (1.25) 

which extends to a bounded bilinear form on H$ x He- L acts by definition isometric from 
D(E) to D(&), with respect to the norms of He and H$, and it extends to an isometric 
isometry from He to H$. Thus is identified with the dual space of He naturally, via the 
extension of the L 2 pairing (■,•). 

Remark 1.1. To give some insight into the interpretation of the spaces He and H§, let us 
mention how He can be characterized in the case of the Newtonian potential |13| : 

H E = {/ : R d -> R measurable : 

there exists an E'-Cauchy sequence f n £ D(E) with f n —tf a.e.}/ ~^ (1.26) 

where ~d is the usual identification of functions that are equal almost everywhere when d > 3. 
For d = 2, in contrast identifies functions that may differ by a constant almost everywhere. 
(It is therefore sometimes said that the massless free field does not exist in two dimensions, 
but that its gradient does. The massless free field is the free field corresponding to <3? in the 
terminology explained below.) To understand this distinction, take a smooth cut-off function 
ipi on M 2 , e.g. with ipi = 1 on -Bi(O) and <p\ = on i?2(0) c , set (p n ( x ) = 9 9 i( :r / n ); an d note 
that E(ip n , (p n ) = n d ~ 2 E(ipi, ipi). Thus, (<p n ) is bounded in He whenever d < 2, and then (by 
the Banach-Alaoglu theorem) there is ip £ He such that cp n — )• ip weakly along a subsequence 
in He] however, <p n — > 1 pointwise, so that ip = 1 G Now -E(l, 1) = implies that the 
constant functions must be in the same equivalence class as the zero function. 

It is well-known that any separable real Hilbert space (H, (•,•)#) defines a Gaussian process 
indexed by H [32]. This is a probability space (fl, P) and a unitary map f £ H (f,4>) £ 
L 2 (P) such that the random variables (/, 0) are Gaussian with variance (f,f )n- Note that 
(/, 0) is merely a symbolic notation for the random variable on L 2 (P) that corresponds to 
f £ H. It cannot in general be interpreted as the pairing of / £ H with a random element 
4>(uj) £ H defined for uj £ £1; see e.g. [30] . 

In particular, if (H,(-,-)jj) is the Hilbert space (i?$, (•, ■)#*)> this process is called the 
/ree yie/d or the Gaussian free field (corresponding to Dirichlet form E or Green's function 
$). The importance of free fields in statistical mechanics, and probability theory in a wider 
sense, is well-recognized. For instance, observables of many models of statistical mechanics 
are intricately related to them, by relations such as the the Kac-Siegert transform [3]. These 
models include spin models such as the Ising model, as well as Coulomb and dipole systems. In 
a different direction, if E is a Markovian form that satisfies some regularity conditions, there 
exists an associated Markov process [19], and it turns out that there are strong connections 
between the distributions of the local times of this Markov process and the free field associated 
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to the same Dirichlet form; see e.g. [Ml [TBI [151 (HJ- 111 particular, in a generalized "non- 
commutative" notion of Gaussian processes that are supersymmetric, this correspondence 
becomes especially striking; see e.g. the review [Tj. The last mentioned correspondence is the 
point of departure for an analysis of the critical behavior of models of self-avoiding walks in 
dimension four [9]. 

For typical applications to statistical mechanics, the measure space (X,(i) of Section [1.31 
is endowed with additional structure such as a distance function, a notion of smoothness, 
etc. as is the case for the Newtonian potential. The global properties of the free field are of 
special interest for statistical mechanics. An example of such a global property is, if X is an 
infinite graph, and X n f X is an increasing sequence of finite graphs approximating X, in an 
appropriate sense, the behavior of 

f Yl e~ v ^ dP(</>), asn^oo (1.27) 
J xex„ 

for some V : R — > R. The covariance $ is typically long-range as in (jl.ip . This makes the 
analysis of the global properties of free fields difficult. 

Decompositions like (|1.16p give rise to notions of scale and corresponding multiscale decom- 
positions of the Gaussian free field and therefore provide a point of departure for multiscale 
analysis. One instance of such an application is the renormalization group method; see e.g. 
[3] and references therein. 



1.5 Main result 

Let (X, fi) be a metric measure space. In addition, let d : X x X — > [0, oo] be an extended 
pseudometric on X. (Extended means that d(x, y) may be infinite and pseudo that d(x, y) = 
for x 7^ y is allowed. Example H31 below gives an example of interest where d is not the metric 
ofX.) 

Let E : D(E) x D(E) — > R be a regular closed symmetric form on L 2 (X) as in Section [1.31 
and denote by L : D(L) — > L 2 (X) the self-adjoint generator of E. Theorem 11.21 assumes that 
(X, fi,d,E) satisfies one of the following two finite propagation speed conditions that we now 
introduce: For 7 > 0, B > 0, and an increasing function 9 : (0, 00) — > (0, 00), let us say that 



(X, fi,d, E) satisfies (P 7) e) respectively (Pq b ) if: 



supp(cos(L2 7 t)u) C iV 6 )( t )(supp(tt)) for all u G C C (X), t > 0, (P 7 



respectively 



E(u,u) < £?||u||£2pn f° r a h u G L 2 (X), 

supp(L n u) C iV e ( n )(supp(u)) for all u G C C (X), n G N, 



e,BJ 



where as before N t (U) = {x G X : d(x, U) < t} for any U C X. The left-hand side of (P 7t $) 
is defined in terms of functional calculus for the self-adjoint operator L. 

Note that if L = — A K d is the Laplace operator of R d , then u(t,x) = [cos(L^t)uo](x) is a 



solution to the standard wave equation (II. 131) . and the condition (Py t e ) with 7 = 1 and 0(t) = t 
is the finite propagation speed property (I1.14p . The property holds for more general elliptic 
operators and elliptic systems (not necessarily of second order), however; see Example 11.41 
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below. Similarly, if L = —A Z d is the discrete Laplace operator (|1.15p . then (Pa B ) holds with 
B = 2 and 9{n) = n, since Lu(x) only depends on u(y) when x and y are nearest neighbors. As 



for the property (-P 7) e), the condition (Pg B ) remains true for more general discrete Dirichlet 



forms; see Examples 1 1 . 4ffL5l 



Let us introduce a further condition: The heat kernel bound (H ayU) ) holds when the heat 
semigroup (e~ tL ) t >o nas continuous kernels p t for all i > and there is a > and a bounded 
function uj : X — > M + such that 



Pt 



(x, x) < uj(x)t- a/2 for all x£X. 



{H 



Criteria for (H au ) are classic; see e.g. [26J for second-order elliptic operators and also the 



discussion in the examples below. 



Theorem 1.2. Suppose (X, (x,d, E) satisfies (P-y$ 
form (jl.24p admits a finite range decomposition 
$t are bounded quadratic forms with 



Mf,9)\ < ^ 2/ l/lbmllslbpo for all f,ge L 2 (X) 



or (Pa B ). Then the corresponding Green 
2]) with S = $ and St = &t such that the 



(1.28) 



Moreover, (H a ^) implies that the <I>t have continuous kernels 4>t that satisfy 



h(x,y)\ < C aa ^uj{x)io(y)t- {a - 2 ^\ 



(1.29) 



1.6 Examples 



Example 1.3 (Elliptic operators with constant coefficients). Let a = (a-ij)i j=i ...,d be a strictly 
positive definite matrix in M. dxd and 



d 

E a (u,v)=y" / {Diuix^aiADMx)) dx, u,v G C™(R d ), 

u,v G C c (Z d ), 



d 



KM= J2 S(Vit*(o;))oi i (V i t;(a;)), 

where Diu(x) is the partial derivative of u(x) in direction i = 1, . . . , d, 

Viu(x) = u(x + ei) — u{x) 



(1.30) 
(1.31) 

(1.32) 



with ei the unit vector in the positive ith direction, and C c (Z d ) is the space of functions 
u : Z d — > R with finite support. For m 2 > 0, further set 



E a ,m?{u,v) = E a (u,v) +m 2 u(x)v(x) dx 



(1.33) 



and define E* m i analogously. Assume that the eigenvalues of a are contained in the inter- 
val [B 2 ,B^], and in the discrete case also that m 2 G [0,M+], for B 2 _,B\,M\ > 0; these 
assumptions are only important for uniformity in the constants below. 

In the continuous context, let d be the Euclidean distance on X = M. d and [i be the 



Lebesgue measure. It follows that (X, n, d, E) satisfies (P 7 g| with 7 = 1, 6(t) = B + t; see 
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Example ll.4l for more details. In the discrete context, let d be the infinity distance on X = Z rf , 
i.e., d(x,y) = maxj = i r ..^ \x{ — yi\, and fi be the counting measure. Then (Pq B ) holds with 
B = B + +Ml and 0{n) = n. 

Theorem 1 1 . 2 1 thus implies that the Green's functions associated to E a ^ m 2 and E* m 2 admit 
finite range decompositions. Let us denote their kernels by (j)t(x, y; a, m 2 ) and 4>l{x, y; a, m 2 ). 
In addition to (jl.29p . it is not difficult to obtain estimates on the decay of the derivatives 
of 4>t and 4>ti like (|1.12p . in this situation of constant coefficients. Since these estimates are 
of interest for applications, we provide the details in Section 13.21 (in a slightly more general 
context). We show that there are constants > depending only on i?_ and B + , and in 
the discrete case also on M + , such that 

\D l a a D 1 ^ D l «D l *ct) t (x,y,a,rn 2 )\ < c^ k t- {d ~ 2) ~ lx ~ ly+2lm2 (1 + m 2 t 2 y k (1.34) 

and 

{D^D^V^V 1 * (j)* t (x,y,t;a,m 2 )\ < Q :k t~ {d ~ 2) ~ l *- l y +2lm2 (1 + m 2 t 2 y k (1.35) 
for all integers l a , l m 2, l x , l y , and k such that 

l m * < \{d + l x + l y ), (1.36) 

and that the following approximation result holds: There is c > such that 

V^V^*(x,y;a,m 2 ) = D l *D l y y<j> t {cx, cy; a, m 2 ) 

+ Oir^-^y- 1 ^ + m 2 t 2 y k ). (1.37) 

This reproduces and generalizes many results of [U Q] . More precisely, we verify that there 
exists a smooth function cj) : R d x [I? 2 . , x [0, oo) — > R supported in |x| < B+ such that 

Mx, V, a, m 2 ) = t^ d ' 2) 4>(^-- a, m 2 t 2 ) (1.38) 

which has the same structure as (jl.lOj) when m 2 = 0; this is scale invariance. Moreover, by 
(|1.37p . the discrete Green's function has a scaling limit and the error is of the order of the 
rescaled lattice spacing 0{t- v ). This result improves |oJ. 

Example 1.4 (Elliptic operators and systems with variable coefficients). Let M £ N and 
dij : R d —7- R AIxM , i, j = 1, . . . , d, be the smooth coefficients of a uniformly elliptic system (or 
in particular, if M = 1, of a uniformly elliptic operator): 

M d 

B 2 _\e,\ 2 <Y,Yj a ih x ^H l j < K\z\ 2 for a11 £ G RdM > x G Rd i ( L39 ) 

k,l=l i,j=l 

with B_,B + > 0. Let us write u = (u 1 , . . . , u M ) G R dM with u' G R d , i = 1, . . . , M. Let 

d „ 

£(u,v)=V / (D i u k (x))a k Ux)(D j u l {x)) dx, u, v e C c °°(R d , R M ) (1.40) 
and analogously in the discrete case (as in (|1.30p . (jl.3ip ). 
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To apply Theorem 11.21 (X, d) is defined by X = R d x {1, . . . , M}, fi is the product of 
the Lebesgue measure on M. d and the counting measure on {1, ... , M}, and the distance is 
given by d((x,i), (y,j)) = d(x,y). In particular, d is only a pseudometric on X. We may use 
the identification of u : M. d — > M A/ and u : X — > R by u(x, i) = w l {x). 

It suffices to verify the condition (Pi^s+t) f° r smooth, compactly supported uq : M d — > W M . 
For such a Uq, set, by using spectral theory for self-adjoint operators: 



u(t) := cos((L + m 2 )2t)u . 



Then, since uq is smooth, u(t, x) : 



is smooth jointly in (t,x), and 



dfu + Lu + m 2 u = 0, d t u(0) = 0, u(0) = u 



n.4i 



(1.42) 



holds in the classical sense. If M = 1, m 2 = 0, and a is the d x d identity matrix, (P\ t t) is the 
finite propagation speed of the wave equation. 

Similarly, in the general situation, the property (P\ t B + t) can be deduced from the finite 
propagation speed of first order hyperbolic systems. This is well-known, but the explicit 
reduction for the case of (|1.42j) with (jl.40p is difficult is to find in the literature. Let us 
therefore sketch how to convert (jl.42p to a hyperbolic system for readers interested in this 
case. For example, one can define v:Kxl d 4 r(<i+2)m by . 



d M 
j=l 1=1 



u 



l,...,d, k = l,...,M, 



Cl-43) 



It follows that v satisfies 
d 



Sd t v + J2 A 3 9 x 3 v + Bv = 0, 



v(0) = (0,(a J Du ) 1 ,...,(a J Du ) d ),u ) 



[1M) 



where S, Aj, B : 



fi{d+2)Mx(d+2)M are ^gfingd as block matrices 



IjV/xM OdMxM OmxM 



OMxdM a 



-1 



OMxdM 



OmxM OdMxM 1-MxM 



Oixi rfx i m 2 
B = I Oixd Odxd Olxd 
-m 2 dx i Oixi. 



1 



MxMi 



:i.45) 



and 



/ o 






















0/ 



IjtfxM, i = 1, . . . ,d. 



(1.46) 



It is immediate that this system is symmetric uniformly hyperbolic, by the symmetry and uni- 
form ellipticity of the matrix a. The property (Pi t B+t) now follows from the finite propagation 
speed of linear hyperbolic systems; see e.g. [3], [25] . 

Nash showed [26J that (Hd u ) holds when M = 1. In |23[ 124] . conditions are given for 
{Hdu) to hold when M > 1. In particular, this includes the constant coefficient case. The 
latter case can be treated by using the Fourier transform; see Section 13.21 
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Example 1.5 (Random walk on graphs). Let (X, E) be a (locally finite) graph, with vertex set 
X and edge set E C P2PO, where X is a countable (or finite) set and P2PO are the subsets 
of X with two elements. Let d : X x X — > [0, 00] be the graph distance on (X, E), i.e., d(x, y) 
is the (unweighted) length of the shortest path from x to y. 

Suppose that edge weights fi xy = ti yx > 0, x, y G X are given. These induce a natural 
measure, also denoted [a, on X by: 

Hx = y~] fixy, K A ) = y^ y ^x for all AO X. (1.47) 

yeX x£A 

The associated Dirichlet form is 

E(u,u) = i ^ Hxy(u(x) - u(y)) 2 for all « G £>(£) = L 2 0u) (1.48) 

and its generator is given by 

Lu{x) = li^ 1 ^xy{u{x) — u{y)) for all finitely supported u : X — > R. (1.49) 

L is called the probabilistic Laplace operator associated to the simple random walk on the 
weighted graph (X, fj,) with transition probabilities \x X yl l^x- Let us remark that a probabilistic 
interpretation (or a maximum principle) does not hold in general for Examples I1.3H1.4I (when 
a is non-diagonal or vector- valued) . 

The Dirichlet form f j 1 . 48 j) is bounded on L 2 (fi) with operator norm 2 so that the property 



(Pq b ) holds with 9{n) = n and B = 2, and Theorem 11.21 is applicable. 

For applications, it is often useful to add a killing rate to the random walk: The proba- 
bilistic Green density with killing rate k G (0, 1) is defined by: 

G K (x,y) = J>>,2/K = (kL + (1 - «))- x ^,y) = {L K )-\x,y) (1.50) 

n>0 

where p n (x,y) is the kernel of the operator P n on L 2 (fi). Note that (|1.50p only converges 
for k = when the random walk is transient, but that L _1 still makes sense as a quadratic 
form on its appropriate domain when the random walk is recurrent, as in (|1.16p . (|1.17p for 
d = 2. Note further that spec(L K ) C [0, 2] for all k G [0, 1], so that Theorem 11.21 is applicable 
uniformly in k G [0,1]. 

Closely related to the killed Green's function G K is the resolvent kernel of L. The resolvent 
of L is defined on L 2 (/i) by G m 2 = (L + m 2 ) -1 for m 2 > 0. It is related to the killed Green's 
density by: 

G K = k- 1 G {1 _ k)/r . (1.51) 

One difference compared with the killed Green's function is that L + m 2 is not bounded 
uniformly in m 2 > 0. To achieve the condition (P2 B ) for fixed B > 0, it is therefore necessary 
to restrict to m 2 < Ml with M\ = B-2. 

Remark 1.6. Other examples which Theorem 11.21 is applicable to include Dirichlet spaces 
that satisfy a Davies-Gaffney estimate [31] such as weighted manifolds and quadratic forms 
corresponding to powers of elliptic operators like A 2 . 
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1.7 Remarks 

Remark 1.7. Theorem 1 1.2 1 also gives the decomposition into sums as in [TOj, [6l [TJ : Suppose that 
the assumptions of Theorem 11.21 are satisfied and, for notational simplicity, that the resulting 
decomposition has a kernel. Then, for any L > 1, 

9(x,y) = ^2Cj(x,y) for all x,y 6 X x X (1.52) 

where the functions Cj : X x X — > [0, oo), jgZ are given by 

/"^ (it 

Cj(x,y) := <j) t {x,y) — forallx,y€X. (1.53) 

They satisfy the following properties: 

Cj is the kernel of a positive semi-definite form, (1-54) 

Cj(x,y)=0 for all x,y £ X with d(x,y) > L j , (1.55) 



and, if (H atUJ ) holds 



\Cj(x,y)\ < c a (x,y) < 



' L -(«-2)(i-i) ( a>2 ) 

L( 2 " a )j (a < 2) (1.56) 

[log(L) (a = 2) 



with c a (x,y) is independent of L. Thus, (Cj)j G z is a finite range decomposition into discrete 
scales of the Green's function Similarly, gradient estimates such as (jl.34p . (|1 .35 j) . (|1.37j) 
in Example 11.31 have obvious discrete versions. 

Remark 1.8. More generally than in Theorem 11.21 we may consider a family of symmetric 
forms, E s , s £ Y, where Y is a domain in a Banach space, with generators L s . Let us assume 
that E s is smooth in s, in the following sense: There exists a projection-valued measure P on 
a measurable space M and a function V : M x Y — > (0, oo), smooth in Y, such that 

F{L S )= I F{\)dP{=! F(V(s,r))dP T . (1.57) 

Jspec(L s ) J M 

An example of this condition is E s (f, f) = E(f, /) + s(f, /) in which case V(s, A) = A + s 
and (L s ) _1 is the resolvent of L; similarly, the killed Green's function of Example 1 1 . 5 1 can be 



expressed in this way. Then the family of kernels 4> s is continuous in s, and if (H a>UJ ) holds 
for s = 0, and V(X,s) > z 2 (s)V(\, 0) + m 2 (s), then 

\$(x,y)\ <C tt)1)1 ^(iM!,)(2( S )r M/7 (l + m)^ (1-58) 
This can be verified by a straightforward adaption of the proof of Theorem 11.21 
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2 Proof of Theorem IT21 
2.1 Spectral decomposition 

The starting point for the proof is the spectral representation of the Green form (11 .24H : 

$(/, /) = / A" 1 d(f, P x f) for all / G £>($), (2.1) 

J spec(L) 

where / G D($) implies that the integral can be restricted to spec(L) \ 0. The main result 
follows by decomposition of the function A -1 : spec(L) \ — > R+. Different decompositions 
are needed under the two conditions (P 7 ,#), (Pg B ). The main idea of the proof is that 
decompositions with good properties exist. The result that we prove after using it to deduce 
Theorem 11.21 is summarized in the following lemma. 



Lemma 2.1 (Spectral decomposition). Suppose that L satisfies (P-yfl) or (Pg B ); in the second 
case, we assume that 7 = 1. Then there exists a smooth family of functions Wt G C° 
t > 0, such that for all A G spec(L) \ 0, t > 0, and all integers I, 



t~W t (X) 



dt 



W t {\) > 0, 

(i + th) l w t (\) <c h 



and that for all u G C C (X), 



supp(W t (L)u) C Afy t) (supp(u)). 

Remark 2.2. More precisely, we will give explicit formulae for Wt that imply 

gr, 



I \ m 



-W t {\) 



(2.2) 

(2.3) 
(2.4) 

(2.5) 
(2.6) 



for all m and /, improving (|2.4p . This improvement is used in Section 
Proof of Theorem [7TH It follows from (f2~2j) that, for any / G £>($), 



*(/,/) 



spec(L) \J 



t-W t (X)j) d(f,P x f) 



(2.7) 



2 



spec(L) 



W t (X)d(f,P x f)\ j 



tf(f,W t (L)f) 



dt 



The exchange of the order of the two integrals in the equation above is justified by non- 
negativity of the integrand, by (|2.3|) . The latter also implies that (f,Wt(L)f) > for all 
/ G L 2 {X). The polarization identity allows to recover &(f,g) for all /, g G D(Q). Finally, 
(|2.5p completes the verification of (|1.8p for <£ t defined by 



Mf,g) = t-(f,w t (L) g ). 



(2.8) 
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It remains to prove that (pa^l) implies f|1.29j) . The semigroup property and the continuity 



of p t imply that p t G C(,(X, L' 2 (X)) with 



\\Pt(x,-)\\ L 2( X ) = Pt(x,y)pt(y,x) dfi(y) =p 2t (x,x), (2.9) 
\\p t (x,-) ~Pt(y, -)IIl 2 (x) =P2t(x,x) +P2t(y,y) - 2p2t{x,y) -)■ as x -> y. (2.10) 



This implies that e _ti : L 2 (X) — >• Cf,(X) is continuous (since e~ tL f(x) = (pt(x, •), /)). Duality 
then implies continuity of e~ tL : Cb(X)* — > L 2 (X) (with respect to the strong topology on 
C{,(X)*). Let M(X) C Cb(X)* be the space of signed finite Radon measures on X equipped 
with the weak-* topology. Let rrii £ M(X) with m; — > 0. Then: 



x 2 \ i 

/ r / p \ \ 2 



|e m\\L 2 (X) = [ l(j x Pt(x,y) drrii(y) } dp,(x) 



(Pt(y, -),Pt(z, •)) dmj(y) dmi(«) -> (2.11) 
ix Jx J 

which means that e~ tL : M(X) — > L 2 (X) is continuous (because X is separable and therefore 
the weak-* topology of M(X) is metrizable). This implies that (l+i 2/7 L)^ : M(X) -> L 2 (X) 
is likewise continuous for all I > a/4. To see this, we use the relation 

(l + t 2 ^X)- 1 =T(l)- 1 / e -V-V s * /7A ds (2.12) 
Jo 

which holds by the change of variables formula and the definition of Euler's gamma function. 
The spectral theorem thus implies that, for any u £ L 2 (X), 

\\{l + t 2 hLT l u\\L* { x)<ni)- 1 / e- s s l - l \\e- st2hL u\\ L , {x) ds. (2.13) 

J o 

Since \i has full support, L 2 (X)nM(X) is dense in M(X) (where L P (X) is always with respect 
to n), and the claimed continuity of (1 + t 2 ^L)~ l : M(X) -> L 2 (X) follows from (T2TTTT) . In 
particular, the pointwise bound for p t implies that for I > a/4, 

/■ oo 

11(1 + ^/7^-/^^ < r(/) -i / e-V-^e-^^H^m ds (2.14) 

JO 

/•oo 

< T(l)- l ^J{x)t-~ a ' 2 ' 1 / e^s'- 1 -"/ 4 ds 







= Cv/^)r Q/27 . 

Let Kt(X) = Wt(X) 1 ^ 2 . Then (12.41) and the spectral theorem also imply that 

\\K t (L)(l + t 2 ^L) l \\ L 2 {x) _ L 2 {x) = sup« t (A)(l + t 2 ^X) 1 < C u (2.15) 

A>0 

uniformly in t > 0. It follows from (f57H|) that K t (L) : M(X) -)• L 2 (X) with 



M^Hl* < C7^)t- a / 27 . (2.16) 
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Finally, by the Cauchy-Schwarz inequality, 

\Mx,y)\ = t 2 ^{K t (L)5 y ,Kt{L)6 x ) < t 2 ^\\ Kt (L)5 y \\ L 2 {x) \\ Kt (L)5 x \\ L 2 {x) (2.17) 

which, with (I2.16p . proves (ll.29p . The continuity of <f>t is implied by the continuity of Kt(L) : 
M(X) -)• L 2 (X) and of 5 X in x £ X (in the weak-* topology). □ 

Remark 2.3. The decay for ffi claimed in (|1.58p can be obtained by a straightforward gener- 
alization of the above argument, replacing (|2.12p by 

(1 + t 2 hz 2 \ + t 2 ^m 2 )- 1 = Til)- 1 / e -s s l-l e -st*hm-> e -szH*h x dg _ (218) 

Jo 

Remark 2.4. Furthermore, by (|2.4p . the operators Wt(L) are smoothing for t > 0, in the 
general sense that, for any t > 0, 

oo 

W t {L) : L 2 (X) -)• C°°(L), where <7°°(L) := p| D(L n ) C L 2 (X) (2.19) 

n=0 

is the set of C°°-vectors for L; see [28]. Standard elliptic regularity estimates imply e.g. that 
C°°{L) = C°°(X) when E is the quadratic form associated to an elliptic operator with smooth 
coefficients. 



2 . 2 Proof of Lemma 1231 

To complete the proof of Theorem II. 2^ it remains to demonstrate Lemma 12.11 We first prove 



it under condition (Pj t g) in Lemma 12.51 below; this proof is quite straightforward using the 



assumption and (jl.2p . Subsequently, we prove Lemma l2~T1 in the situation of condition {Pq B ) 
in Lemma \2.7\ here additional ideas are required. 

To fix conventions, let us define the Fourier transform of an integrable function <p : R — >■ R 

by 

0(k) = (2TT)- 1 / tp(x)e- ikx dx for all fcet (2.20) 

Lemma 2.5 (Lemma 12.11 under (P~g|). For any ip : R — > [0, oo) such that (p is smooth and 
symmetric with supp((^) C [—1,1], and for any 7 > 0, there is C > such that 



W t (X) := Cip{\^t) 



(2.21) 



satisfies (p!2"]) . (|231) . P^l . and a/so ([12]), /or all \>0,t> 0; and i/ ( |7^| ) /ioZds, i/ten (W t ) 
a/so satisfies (|2.5[) . 



Remark 2.6. It is not difficult to see that such exist. For example, if k is a smooth real-valued 
function with support in then <p = \k\ 2 satisfies the assumptions. For simplicity, let 

us assume sometimes in the following that ip is chosen such that C = 1 when Lemma 12.11 is 
applied. 
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Proof. Note that for any ip : [0, oo) —> M with t<p(t) integrable, there is C > such that 



\- 1 = C / t^(p(X^t) — for all A > 0. (2.22) 
Jo t 

This simply follows (as in (jl.2p ) because the right-hand side is homogeneous in A of degree 
— 1, which is immediate by rescaling of the integration variable. This shows (|2.2p ; (|2.3p is 
obvious by assumption; and (12. 4p follows since (p is smooth. The improved estimate (|2.6I) 
follows from the chain rule (or Faa di Bruno's formula) and 



d m i 
■A2 7 



^ Cy,m (2.23) 



for non- negative integers m, using that supp(<^) C [—1,1] implies that ip is smooth. Moreover, 
since supp((/i>) C [—1, 1], and since (p is smooth, 

W t (L)u = C J 0(s) cos(L^ts)u ds for all u G L 2 (X), (2.24) 

where the integral is the Riemann integral, i.e., the strong limit of its Riemann sums (with 
values in L 2 ). Therefore (|2.5p follows from (Pyg). □ 



The previous proof makes essential use of the finite propagation speed of the wave equation 
(P~/fi) to prove (|2.5p . This property fails for discrete Dirichlet forms such as (jl.3ip where we 



instead know the property (Pag) that polynomials of degree n of the generator have finite 
range 9{n). 

This leads to the following problem. Find polynomials Wj*, t > 0, of degree at most t 
satisfying the properties (|2.3p . ([2.4]) . (]2.6j) such that the decomposition formula ()2.2[) for 1/A 
holds. In the proof of Lemma l2.5] the verification of (12. 4p (and (|2.6p ) and of the decomposition 
formula (]2.2p are directly linked to the "ballistic" scaling of the wave equation: Wt(X) = 
VF^At 2 ). To construct polynomials satisfying such "ballistic" estimates, we are led by the 
following remarkable discovery of Carne [11]: The Chebyshev polynomials T^, k G Z, defined 
by 

T k (6) = cos(A;arccos(<9)) for all 6 G [-1, 1], k € Z, (2.25) 

are solutions to the discrete (in space and time) wave equation in the following sense: Let 
V+/(n) = f(n + 1) — f(n) and V_/(n) = f(n — 1) — f(n) be the discrete (forward and 
backward) time differences. Then, as polynomials in X, 

V_V + T n (X) = V+V_T n (X) = 2{X - l)T n (X). (2.26) 

In particular, when2(X — 1) = —L or equivalent ly X = 1 — thenf(n,x) = [T n (l — ^L)u](x) 
solves the following "Cauchy problem" for the discrete wave equation: 

-V + V-v + Lv = 0, v(0) = u, (V+u- V_u)(0) = 0. (2.27) 

The analogy between the discrete- and the continuous-time wave equations is like that be- 
tween the discrete- and the continuous-time random walk. It turns out that the structure of 
Chebyshev polynomials allows to prove the following lemma. 
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Lemma 2.7 (Lemma 12.11 under (Pa B ))• Let ip : M — > [0,oo) satisfy the assumptions of 
Lemma POl Then W£ : [0,4] — > [0, oo), defined by 

W t *(X) := J^(arccos(l- \\)t-2imt) for all A G [0,4], t > 0, (2.28) 

is the restriction of a polynomial in A of degree at most t to [0,4], with coefficients smooth in 
t, and, for any e > 0, $M$), $ZM, and $2J^ hold for all A G (0,4 - e], t > 0. 

Proof. The proof verifies that W t * as defined in (|2.28[) has the asserted properties. Let 

tp*(x) := ^2 ~ 27rnt ) = ^2 t ~ 1( f( k / t ) cos(kx) (2.29) 

n&L k€Z 

where the second equality follows by symmetry of <p, the change of variables formula, and 
a version of the Poisson summation formula which is easily verified, for sufficiently nice (p. 
Then the claim (12, 2p can be expressed as 

poo J± 

A" 1 = / iV£(arccos(l - \\)) — for all A G (0,4]. (2.30) 
Jo " t 

Let x = arccos(l — ^A) or equivalently A = 2(1 — cosx) = 4sin 2 (ix). In terms of this 
change of variables, (|2.30p and thus the claim ()2.28j) are then equivalent to 

r°° dt 

isin- 2 (ix) = / t 2 if* t (x) — for all x G (0, it}. (2.31) 
Jo t 

The left-hand side defines a meromorphic function on C with poles at 27rZ. Its development 
into partial fractions is (see e.g. (2J page 204]) 



\ sin" 2 (ix) = - 27m) -2 for all x G C \ 2vrZ. (2.32) 

It follows, by (|2.22p with 7 = 1 and A = (x — 27m) 2 , assuming C = 1, that 

poo ^ 

±sin _2 (±a;) = V / t 2 tp((x - 2vrra)t) — for all x G (0,tt]. 



(2.33) 



The order of the sum and the integral can be exchanged, by non-negativity of the integrand, 
thus showing (|2.3ip and therefore (|2.2p . 

To verify that W± is the restriction of a polynomial, we note that by (|2.28j) . (I2.29p . and 
supp((£) C [-1,1], 

W t *(X) = c^(arccos(l - ±A)) = ^ t~ 1 fi(k/t) cos(k arccos(l - |A)) (2.34) 

= Yl t- l £(k/t)T k (l - ±\) 
keZn[-t,t] 

where T^,/c G Z, are the Chebyshev polynomials defined by (|2.25|) . This shows that W*(X) 
is indeed the restriction of a polynomial in A of degree at most t to the interval A G [0,4]. 
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In particular, (|2.5p is a trivial consequence of (P$ B ) which states that polynomials in L of 
degree n have range at most 9{n). 

Finally, we verify the estimate (|2,6p and thus in particular (|2.4p . To this end, we note 
that, in analogy to (|2.23[) . for AG [0,4 — e] and non-negative integers m, 



A m "2 



arccos(l — ^A) 



<C £ 



(2.35) 



For example, for m = 1, 



— arccos(l - ±A) = ±(A- ^A 2 )^ < £ -|a^ for A € [0,4 - e]. (2.36) 

OA 



2'V ~~ 2 

Therefore (|2.6p follows, by the chain rule (or Faa di Bruno's formula), from 

gm 



(1 +r(l -cos(x))'t 



0X r 



< C, 



1.771 



(2.37) 



which we will now show. The argument is essentially a discrete version of the classic fact 
that the Fourier transform acts continuously on the Schwartz space of smooth and rapidly 
decaying functions on R. To show (|2.37p . first note that 



(1 - cos(x))e ikx = e ikx - l e i( - k+1)x - \e i{k ~ l)x =: A k e ikx 
and thus by induction, for any I E N, 

(1 - cos(x)) l e ikx = (1 - costs)) 1 - 1 Afce** = A fc (l - cos(x)) 1 - 1 e lkx = A l k e ikx . 
It follows by (12.29P and summation by parts that 

gm 



{l + t 2 (l-cos(x)yt 



lj.—m 



dx^^ (X 



) = ^t- 1 <p(k/t)(ik/t) m [(l + t 2 A k ) l e 



ikx] 



(2.38) 



(2.39) 



(2.40) 



k&Z 



ikx 



k& 



Let h(s) = T%(\s\ — l)l| s |<i for s£R. Then, for any smooth / : 

A n k f(k) = (h* n *D 2n f)(k), 



(2.41) 



where * denotes convolution of two functions on R, h* n = h * h * ■ ■ ■ * h, and Df is the 
derivative of /. Indeed, 

A fc /(fc) = ~\ f [Df(k + t) - Df(k - t)} dt (2.42) 
Jo 

= -\l [ D 2 f(k + s)dsdt= [ D 2 f(s)h(s-k)ds = (h*D 2 f)(k), 

Jo J-t Jr 

and (|2.4ip then follows by induction: 

A£ +1 / = A{h* n * D 2n f) = h* D 2 {h* n * D 2n f) = h*h* n * D 2 D 2n f. (2.43) 
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It then follows using the facts that ^2 keZ \ h* n (k — s)\ < C n , uniformly in s G R, and that (p 
is smooth and of rapid decay, 

t-^la + t^lYimWk/ty^ (2.44) 

= E C ^E / \h* n (k-s)\\[D 2n ((-r0)](s/t)\ds 
n=o kez jR 

<E C ^ / \[D 2n ((-rms/t)\ds 

n=0 ^ 

/ \[D 2n ((-) m 0)]( s )\ ds<c mtl 

n=0 ^ 



and thus ()2.37|) . and therefore (j2.6|) . follow from this inequality and ()2.40p . □ 



Proof of Lemma \2. 1\ Lemma 12.11 under (P^.g) follows immediately from Lemma 12.51 under 
<P7b1>) [t follows from Lemma O by setting W t (A) = W t *(^A). □ 



3 Extensions 

3.1 Discrete approximation 

In view of the discussion about Chebyshev polynomials before Lemma [2. 71 it is not surprising 
that the functions W£ of Lemma 12.71 approximate the Wt of Lemma 12.51 In Proposition 13.11 
below, we show that this is indeed the case with natural error 0{t~ l ) as t — > oo. This result 
is used in Section 13.21 to prove (I1.37P . 



Proposition 3.1 (Discrete approximation). Let <p be as in Lemma [2.5\ and 2.1 , with associated 
functions Wt and W£ for 7 = 1. Then, for any integer I, 

\W t *(X)-W t (X)\<C l (l\/t)- 1 (l + t 2 X)- 1 for all X € [0,4]. (3.1) 
In particular, Wf(X/t 2 ) — > Cip(Xz) as t — ^ oo. 



Proof. Note that it suffices to restrict to t > 1, since for t < 1, the claim follows from (|2.4p . 
The left-hand side of f|3. 1 j) is then proportional to the absolute value of 

i^(arccos(l - ±A)t) - ip(xh) + V] 9?(arccos(l - §A)i + 2nnt). (3.2) 

nSZ\{0} 

We estimate the difference of the first two terms in (|3,2p and the sum separately, and show 
that each of them satisfies (13.11). The first two terms can be written as 



<^(arccos(l - ±A)i) - (p(xh) = (arccos(l - |A) - X^)t( t (X) (3.3) 

with ^ 

Ct(A) = / <p'(s arccos(l - \X)t + (1 - s)Xn) ds. (3.4) 
Jo 
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The bounds 

V2X = arccos(l - A) + O(A) as A ^ 0+, (3.5) 

V2\ < arccos(l - A) < f V2A for all A G [0, 2], (3.6) 

and the rapid decay of ip' therefore imply that 

|Ct(A)| <Q(1 + Xt 2 )- 1 (3.7) 

and 

<£>(arccos(l - ±A)i) - ip(\h) < Cjt -1 (l + i 2 A)^. (3.8) 

To estimate the sum in (|3.2h , we can use the rapid decay of (p with the inequality x + y > 
2(xy) 1 / 2 to obtain that 

E V?(xt + 2vrnt) < C x ^ (1 + xt + 2-Knt)~ l (3.9) 
nez\{o} neZ\{o} 

< Q(l + xt)-'/ 2 t-'/ 2 £ n^/ 2 < £7,(1 + xt)-'/ V'/ 2 

n>0 

for any / > 2, with the constant changing from line to line. In particular, upon substituting 
x = arccos(l — ^A), this bound and (|3.6|) imply 

^(arccos(l - ±A)t + 2imt) < C^~ 2i (l + t 2 X)' 1 . (3.10) 

ngZ\{0} 

The claim then follows by adding (|3.8p and f)3. lOj) . □ 



3.2 Estimates for systems with constant coefficients 

In this section, we verify the assertions of Example 11.31 We work in the slightly more general 
context of second-order elliptic systems (instead of operators) with constant coefficients. These 
are defined as in Example 11.41 and the claims of Example 11.31 hold mutadis mutandis. The 
analysis is straightforward, with aid of the Fourier transform. It reproduces results of [I]. 



3.2.1 Spectral measures 

The spectral measures corresponding to the vector-valued case of (|1.30j) are given in terms of 
the Fourier transform as follows. For F : [0, oo) — > R, 



M .. 

(v,F(L)u)= £ / 
k,l=i Jr 



(3.11) 



ki 



where u = (u 1 , . . . , u ) is the component-wise Fourier transform of u = (it 1 , . . . , u M ), 



(3.12) 



kl=l M 
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are symmetric positive definite M x M matrices, for all £ E R rf , and the matrices F(a(£)) are 
defined in terms of the spectral decomposition of a(£). Similarly, for the (vector- valued case 
of the) discrete Dirichlet form (|1.31|) , 



M . 

(v,F(L)u)= / 

k,i=i J l-*>*. 



1 


H 


>[-TT,1T] d 





*>W(0 dt (3.13) 



where here u is the component-wise discrete Fourier transform. Let us also write 



«*(£) := E "v n - e ^)(! - c_<& ) = E 4K 1 - - < 



v»J=l 



feJ=l,...,M 



(3.14) 



We will often use, without mentioning this further, that the spectra of a(£) and a*(£) are 
bounded from above and from below by |£| 2 . 



3.2.2 Estimates 

Let us introduce the following notation for derivatives: For a function u : R d — > R, we regard 
the Zth derivative, D l u(x), as an Z-linear form, and \D u(x)\ is a norm of the form D l u{x). In 
terms of the Fourier transform, we denote by the corresponding "multiplier" operator 

from functions to /-linear forms, and by |l)'(£)| its norm. Similarly, for a discrete function 
n : Z rf — t- R, the Zth order discrete difference in positive coordinate direction is denoted by 
V'tt(x) and has Fourier multiplier V'(£). In particular, when I = 1, 

V(£)-(e^-l,...,e^-l). (3.15) 

Furthermore, and p will denote integers that may be chosen arbitrarily, and C constants that 
can change from instance to instance and may depend on k and p, as well as Z = (l x , l y , l a , l m 2 ), 
B + , and M + , but not on x, £, and m. 

Proof of (|1.38p . (ll.34p . ([1.35p . It follows by the change of variables from the fact that 

a(£) is homogeneous of degree 2, and from Wf(A) = W\(\t 2 ) that 

y; a, m 2 ) = t 2 [ W t (a(0 + m 2 )e i(x-,,H d £ ( 3 - 16 ) 
= t-("- 2 )^(^;a J m 2 t 2 ) 

with 

#r;a,m 2 ):= / Wi(a(£) + mV * - ^ d£ (3.17) 

which is supported in |x| < B + . This verifies (|1.38p . Furthermore, (|1.34p is a straightforward 
consequence of f|3. 16j) by differentiation and (|2.6p . Let us omit the details and only verify 
them explicitly in the discrete case (|1 .35[) : The (derivatives of the) decomposition kernel §\ 
can here be expressed as 

i#2?& 3 V|rVj»#fo^ (3.18) 
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with 

<^(x;a,m 2 ) = t <*+i*+i w -2f m 2 /" D l a a Dfy 2 W? (a* (£) + m 2 )^V^e^^ df. (3.19) 

•/ [— 7T,7r] d 

Thus (H2D, |V(OI < C\£\, and • a*(£)»7 > C||| 2 |?7| 2 for 77 € M M imply 

|<^(x;a,m 2 )| < C / (1 + C|^| 2 t 2 + m 2 t 2 )- k - p (t\Z\) l * +l v- 21 ™ 2 t d d£ (3.20) 

V [-7T,7r] d 

< C(l + m 2 t 2 )- fe / (1 + c|£| 2 )-P|^ + ^ m2 

and therefore that the integral converges if ^(d + l x + Z^) > Z m 2 and p is chosen sufficiently 
large. It follows that 

|^(:r;a,m 2 )| < C(l + m 2 t 2 )" fc (3.21) 

verifying the claim. □ 

Proof of (jl.37p . Let us assume that 5 = 3. Then 

V£vi»#(aM/) - L%Dj>Mx,v) = t 2 I W t *(a*(0)^ ly e iHx ~ y) d£ (3.22) 

J{-Tr,w] d 

-t 2 f W t (a(0)D lx & y e iHx ~ y) d£. 

To simplify notation, we will write & = D lx D ly = D lx <g) D ly if I = (l x , l y ), and similarly for 
V. Then the difference (|3.22p may be estimated as follows. Proposition 13.11 implies 

/ \W t *(a*(0 + m 2 ) - W t (a*(0 + m 2 )\\D l (0\ 

J[-7T,7r] d 

<cr l f {i + c\z\h 2 + mh 2 )- p - k \t\ l dt<ct- d - l -\i + mH 2 )- k (3.23) 

where we have assumed in the second inequality above that p was chosen sufficiently large so 
that the integral is convergent. Similarly, we may proceed for the other differences, always 
choosing p large enough in the estimates. Using ()2.6[) with m = 1 and |a*(£) — = 0(|f | 3 ), 
which follows from Taylor's theorem, we obtain 

/ \W t (a*(0 + m 2 ) - W t (a(0+m 2 )\\D l (0\ d£ 

J[-iT,ir} d 

<C I If 1(1 + C\£\ 2 t 2 + m 2 t 2 )- p - k \£\Ui<Ct~ d - l ~ l {l + m 2 t 2 )- k . (3.24) 
Taylor's theorem similarly implies that |V*(f) - < C|f | m so that, by ([211) . 

[ \ Wt *(a*(0+rn 2 )\\V l (0-D l (0\dt 

J{-n,w] d 

<C I {1 + C\£\ 2 t 2 + m 2 t)- p - k \£\ l+1 d£ < Ct~ d ~ l ~ 1 (l + m 2 t 2 y k . (3.25) 

22 



Finally, we obtain by (|2.4p that 




W t (a(0 + m 2 )\\D l (0\dd 



< C 




(1 + C|^| 2 t 2 + mH 2 y p - k \£\ l d£ < Ct~ 2 P(l + m 2 t 2 )- k . (3.26) 



The combination of the previous four inequalities gives (jl.37p . 



□ 
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